a: 

Oh! 

43 



CERN-TH/99-301 

HEPHY-PUB 722/99 

UB-ECM-PF 99/15 

UWThPh-1999-60 



The Heavy Quarkonium Spectrum at Order 

mal In a s 

On, 

Nora Brambilla 1 , Antonio Pineda 2 , Joan Soto 3 and Antonio Vairo 4 

O i ' lmti " f ™' che *** Bol tZm ann S as,e 5, A-,090 Vienna, Austria 

^" ■ Theory Division CERN, 1211 Geneva 23, Switzerland 

q^ ' 3 Dept. d'Estructura i Constituents de la Materia and IFAE, U. Barcelona 

Diagonal 647, E-08028 Barcelona, Catalonia, Spain 

Institut fur Hochenergiephysik, Osterr. Akademie d. Wissenschaften 
Nikolsdorfergasse 18, A-1050 Vienna, Austria 

en 

O 



PACS numbers: 12.38.Bx, 12.39.Hg, 14.40. Gx 



Abstract 



We compute the complete leading-log terms of the next-to-next-to-next-to-leading- 
order corrections to potential NRQCD. As a by-product we obtain the leading 
logs at O{mo^) in the heavy quarkonium spectrum. These leading logs, when 
Aqcd <C ma>1, give the complete 0(ma^ In a s ) corrections to the heavy quarko- 
j^ ■ mum spectrum. 
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1 Introduction 



The last years have witnessed important progress in the theoretical understanding of heavy- 
quark-ant iquark systems near threshold through the use of effective field theories [1-7] (see 
also [8-11] for related work in QED). The key point relies on the fact that, since the quark 
velocity v is a small quantity, v -C 1, a hierarchy of widely separated scales, m, mv, mv 2 ..., 
is produced in these systems, where m is the heavy quark mass (hard scale). We will call 
mv and mv 2 the soft and the ultrasoft (US) scale respectively One can take advantage of 
this hierarchy by systematically integrating out the scales above the energies one wants to 
be described by the effective theory. 

After integrating out the hard scale, non-relativistic QCD (NRQCD) is obtained [1]. The 
Lagrangian of NRQCD can be organized in powers of 1/m. The matching coefficients of 
NRQCD are (non-analytic) functions of m. After integrating out the soft scale (proportional 
to the inverse of the size of the bound state) in NRQCD, potential NRQCD (pNRQCD) is 
obtained [3]. The Lagrangian of pNRQCD is organized in powers of 1/m and of the relative 
coordinate r (multipole expansion). The matching coefficients of pNRQCD are (non-analytic) 
functions of r. 

The integration of the degrees of freedom is done through a matching procedure (see [2,12,10,7] 
for details). The matching from QCD to NRQCD can always be done perturbatively since, 
by definition of heavy quark, m ^> Aqcd [2,12]. The matching from NRQCD to pNRQCD 
can only be carried out perturbatively when mv ^> Aq C d- We will assume this to be so 
through this work. Therefore, the matching coefficients both in NRQCD and in pNRQCD 
can be computed order by order in a s . The non-analytic behaviour in 1/m appears through 
logs in the matching coefficients of the NRQCD Lagrangian: 

(Tfl \ 

In h B , 
Vh J 

where fih denotes the matching scale between QCD and NRQCD. A typical matching coef- 
ficient of pNRQCD has the following structure^: 

(Tfl \ 

A' In — + A' In fi h r + B' In fir + C ) , 
Vh J 

where fi denotes the matching scale between NRQCD and pNRQCD. V denotes a function 
of r, p and the spin of the particles, which is analytic in the two last operators but typically 
contains non-analyticities in r. We see that the dependence on fi^ cancels here and that two, 
potentially large, logs appear in the pNRQCD matching coefficients: In mr and In fir. They 



1 We have assumed, since it will be so in the present work, that all the fih dependence cancels in 
the matching procedure. This is in general not true at higher orders due to the fact that there are 
terms in the quantum-mechanical perturbation series that also exhibit an explicit dependence on 
the cut-off /i/j. 



represent the leading-log corrections to the pNRQCD Lagrangian. In practice, the leading- 
log dependence on In mr of the pNRQCD matching coefficients can be trivially extracted 
from the leading-log behaviour of the matching between QCD and NRQCD at one loop, 
which is known, see [2,12]. In order to obtain the leading-log dependence on ln/zr of the 
pNRQCD matching coefficients, it is just sufficient to know the behaviour of pNRQCD in the 
ultraviolet (UV) at next-to-leading order in the multipole expansion. Moreover, everything 
can be worked out at some definite order in 1/m. This will be discussed in greater detail 
below, where we will obtain the leading-log dependence of the pNRQCD matching coefficients 
up to 0(l/m 2 ) and at leading order in the multipole expansion. The leading-log running of 
the matching coefficients at O(l/m ) has already been calculated in [6]. 

If we assume Aqcd < mv 2 the leading-order solution to the heavy quarkonium spectrum 
corresponds to a Coulomb-type bound state. In a Coulomb-type bound state, v ~ a s and all 
the terms of the type a s /v have to be resummed exactly. The next-to-leading order (NLO) 
terms are corrections of 0(v,a 8 ), the next-to-next-to-leading order (NNLO) terms are cor- 
rections of 0(v 2 , va s , a 2 ), and so on. These corrections may get multiplied by parametrically 
large logs. The matching sketched above and exploited in the next sections is sufficient to 
obtain the leading-log corrections of the NNNLO terms in the pNRQCD Lagrangian (there 
are not parametrically large logs at NLO and NNLO). Since results are now available for 
some observables up to NNLO order [13], the leading-log corrections of the NNNLO terms 
in the pNRQCD Lagrangian are the first step towards the complete evaluation of the heavy 
quarkonium spectrum at NNNLO order. It is one of the aims of the present work to set up 
the framework for an eventual full NNNLO calculation. For what concerns the relation be- 
tween the heavy quarkonium mass and the pole mass our calculation will allow us to obtain 
the complete 0(mal lna s ) corrections in the situation mv 2 ^> Aqcd- I n the more general 

situation, Aqcd < mv 2 , we will still be able to obtain all the ma^ In and ma^ In - 

m pi 

terms, where the \x dependence of the latter cancels against US contributions that now should 

be evaluated non-perturbatively. These US contributions have recently been computed per- 

turbatively in the MS scheme in Ref. [14] where the \x dependence can already be read 

of. 

We distribute the paper as follows. In sections 2 and 3 we fix the effective field theory 
framework we are working in. In particular we write the pNRQCD Lagrangian in the equal 
mass case up to order 1/m 2 . In section 4 we perform the matching. In section 5 we apply 
the result to the quarkonium spectrum by computing the leading-log corrections at 0(ma^). 
The last section is devoted to some comments and conclusions. 



2 NRQCD 



After integrating out the hard scale m, one obtains NRQCD [1]. Let us, first, write the most 
general Lagrangian (up to field redefinitions) up to 0(1 /m 2 ): 



£nr«od = *{*7°B„ + £- + JL + Cf9 5_1 (1) 

7°(D-E-E-D) 7°S-(D x E-Ex D)l 
+c D g— — 2 + icsg- > — — '- * 

_ —C C^ u _|_ /^ n2/i/«' _i_ 3 r /^ia ^6 /oc 

, ( - T /i^ ( - r ~t~ ^Kj^ v LJ Ur -f- 2 & J a b c IJ.V fia^* va 

+^Wx + ^ W*X + ^tT-^txa + ^TV^TV*. 
m 2 m 2 m 2 m 2 

We have here ^ = ip + x, where ip is the Pauli spinor field that annihilates the fermion and 
X is the Pauli spinor field that creates the antifermion, S = , iD° = ido — gA° and 

?D = iV + gA. We have also included the V 4 /m 3 term in Eq. (1), since it will be, once 
the power counting is established, as important as other 0(l/m 2 ) operators in the NRQCD 
Lagrangian. This Lagrangian is sufficient to obtain the pNRQCD Lagrangian up to, and 
including, the leading-log terms at NNNLO. One can see that this is so by following an 
argument analogous to Rcf. [10]. The coefficients cp, cd, c$, d<i and d% can be found in Ref. 
[2] and dij (i,j = s,v) in [12]. Since d<i and d% do not contain \nm/nh terms, they may be 
neglected for the present analysis. 



3 pNRQCD 



Integrating out the soft scale in (1) produces pNRQCD [3,7]. The pNRQCD Lagrangian 
reads as follows: 



r)-^--^r + ...\S 



m m 2 



£pNRQCD = Tr|st^ -^ + £^-K (0) (^ 
+0 1 " (iD -^ -V}°\r) + . .) o| 
+gV A (r)Ti {O f r ■ES + S t r-EO} + g^^-t-Tr {O f r • E O + O f Or ■ E} 



1 

-1 
4 



--G%G^\ (2) 



where we have explicitly written only the terms relevant to the analysis of the leading-log 
corrections at the NNNLO; S and O are the singlet and octet field respectively. All the gauge 
fields in Eq. (2) are functions of the centre-of-mass coordinate and the time t only. For a 
more extensive discussion we refer the reader to Ref. [7]. 

The functions V are the matching coefficients of pNRQCD. They typically have a non- 
analytic dependence on r. Although in Eq. (2) this dependence is not shown explicitly, 



the matching coefficients V also depend on other parameters such as the mass (through 
logarithms), or operators such as the momentum and the spin (analytically). For the present 
purposes we can approximate Va, Vb and V^ to their leading-order value in a s (although 
the leading- log contributions are known [7]): 

V A = V B = l- V?\r)=(^-C F )^. (3) 

Let us now display the structure of the matching potentials VJ°\ V} 1 ' and V^ 2 **, which are 
the relevant ones to our analysis. 

1) Order l/m°. From dimensional analysis, VJ°\r) can only have the following structure: 

K (o )(r) = _ Cf ^M. (4 ) 



2) Order 1/m. From dimensional analysis and time reversal V^\ can only have the following 
structure: 



VP CfCaDW 



m 2mr 2 

3) Order 1/m 2 . At the accuracy we aim, Vj- 2 ** has the structure 



VP C F D® (1 2 \ C F Dfl 1 Trgggg (3) 



(5) 



m 2 2m 2 I r ' r J 2m 2 r 3 m 2 

+ !^lW> (r) + ^I L • S + ^fl Sl2 ( f ), (6) 

6m 2m r A 4m r J 

where ^12 (r) = Sf-o - ! f cr 2 — cr!cr 2 and S = (Ti/2 + cr2/2. Note that p appears analytically in 
the matching potentials. The power on n to which p appears in the potentials is constrained 
by the power in 1/m. At order 1/m 2 , p can only appear at most at the square poweiQ. 

The matching coefficients, ay s , D s , contain some lnr dependence once higher order correc- 
tions to their leading (non- vanishing) values are taken into account. In particular, we will 
write expressions like <5® (r) lnr. This is not a well-defined distribution and should be under- 
stood as — reg— , which is the Fourier transform of Inl/k (see [10] and references therein). 

Nevertheless, in order to use the same notation for all the matching coefficients, and since it 
will be sufficient for the purposes of this paper, where we are only interested in the depen- 
dence on the logs of the matching coefficients, we will use the expression 5^(r) lnr, although 
it should always be understood in the sense given above. 



2 This result follows from the knowledge of the operators present in the 1/m expansion in the 
underlying theory (NRQCD) and from the use of time reversal. 



Finally we note that the above representation of the potential can be related to others found 
in the literature [15] by the use of the equations of motion and of the identity 



{;.P 2 } + ^ + 4^'(r) 



P 



r r- (r-pjp 



The ambiguity in the form of the potential is due to the freedom to perform time-independent 
unitary field redefinitions, which do not change the spectrum of the theory. One could try to 
fix this freedom by working with a minimal set of independent potentials (just in the same 
way as the matching coefficients of NRQCD are unambiguous for a given renormalization 
scheme once the Lagrangian is written in the minimal form, namely with no higher time 
derivatives). In principle, given a potential, one could obtain its minimal form by using field 
redefinitions. This point will be elaborated elsewhere. 



4 pNRQCD matching coefficients 



In order to obtain the explicit expressions for the different matching coefficients in pNRQCD 
(i.e. ay 3 , D^\ Df\ ...), one has to perform the matching between NRQCD and pNRQCD. 
A detailed description of the procedure can be found in [3,7,10]. Let us first write our results 
for the matching coefficients up to 0(l/m 2 ). The matching coefficients read as follows: 



&Vs =a s (r)\l + (ai + 2^ E f3 ) 



« s (r) 

47T 
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TV 
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D { H = a s (r), 



D 
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d.s 
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D% 



D 
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(2) . 
LS,s " 

(2) 

Sl2,S 



a s (r) |1 - 
a B (r) |l - 



a s nCp_ 17 V A 
3 3~ 

7C A a s 

4 -k 

2C A a s 



In mr 



In mr 



3 7T 

a s (r) <! 1 — Ca— In mr 

IT 



In mr H Cp m/xr > 

3 7r V 2 / J 



(7) 

(8) 

(9) 
(10) 

(11) 

(12) 

(13) 
(14) 



To obtain ay s with the above accuracy (Eq. (7)), it is necessary to perform the matching 
between NRQCD and pNRQCD (at 0(l/m°)) exactly at the two-loop level and with leading- 
log accuracy at the three-loop level, i.e. to compute the static potential to this order. The 



one-loop result was obtained in Ref. [16], the two-loop one in Ref. [17] and the three-loop 
leading log in [6] . f3 n are the coefficients of the beta function and the values of a\ and a<i can 
be found in Ref. [17] (see also [17] for notation). 

For D^ (Eq. (8)) we need to perform the matching exactly at the one- loop level and to 
obtain the leading-log at the two-loop level. The one-loop matching has been done first (to 
our knowledge) in [18] (note that D^ has no tree-level contributions). It gives 



5D[V > 



<X r 



The two-loop leading-log term is obtained from the UV behaviour of a next-to-leading order 
calculation in the multipole expansion to 0(l/m) in pNRQCD (see Fig. 1). 

For the different D^ terms we only need to compute the matching at tree level (along the 
same lines as in Ref. [10]) and the leading- log at one loop. The leading dependence on In mr 
is obtained by taking into account the NRQCD matching coefficients in the vertices when 
matching to pNRQCD at tree level and by setting \Xh = 1/r. The relevant contributions read 
as follows: 

5D^ = a s (r), 
6D$> a =a B (r), 



SD (2) s = a s (r)(2 + c D -2c 2 F ) + ' 



(LvS \ J^Stl T _, [flsS \ ids 



f a a (2C F 17C A \ 
~a s (r) H I — 1 — J In mr 

SD$ = a s (r)c 2 F - (d sv + C F d vv ) ~ a s (r) \ 1 -^— In mr L 

LT\L> F l- 4 71 J 

xn ( 2 ) a B (r) f 2C A a s \ 

$ d ls,s = ~^-( c s + 2c f) - «s(r) |1 — — In mrj , 

5D { sl s = a s (r)c 2 F ~ a B (r) I 1 - C A — In mr I . 

The leading dependence on In fir is obtained from the UV behaviour of a next-to-leading 
order calculation in the multipole expansion at 0(l/m 2 ) in pNRQCD (see Fig. 1). It is worth 
noting that the spin-dependent matching potentials do not have ln/zr contributions. This 
follows from the fact that the vertices in Fig. 1 have no spin structure. In fact, they are 
known with better accuracy than the one demanded here (see [5,15] and references therein). 

In order to obtain the dependence on fj, of the above matching coefficients, it is sufficient, 
to compute the UV divergences of pNRQCD at the next-to-leading order in the multipole 
expansion and up to 0(l/m 2 ). Basically, we have to compute the diagram shown in Fig. 1 
(see [7] for the Feynman rules in pNRQCD), say 

i oC F C A i f d 3 k i , % 

-r I - — —k 



E-p 2 /m-Vs {0) 3 ' J (2tt) 3 e - k - U o (0) - p 2 /m E - p 2 /m - U s (0) ' 



£ 




%> 



1/(E- 7j°) - p 2 /m) 

Fig. 1. The UV divergences of this diagram in pNRQCD fix the jjl dependence of ay 3 , D s , D s , 
Z?\zP a ndzP. 



where all the derivatives are assumed to act on the right. The divergences (the fi dependence) 
we find in this expression must cancel against the \x dependence of the matching potentials 
in the expression: 

_^l/2 I Z l/2 f 

s E-Vl 0) -p 2 /m-Vs {1) /m-Vl 2) /m 2 * ' 

fixing the fi dependence of the different matching potentials listed in Eqs. (7)-(14). We also 
obtain the leading In fir dependence of the singlet normalization factor Z X J 2 defined in [7], 



Z s,US - V 7Vc {*■ + Z *,US + ~ Z s,XJS + ^ Z S,USP J » {^) 

Z% = ^C F C\al\n^ (16) 

Z% = I i^C 2 F + 2C f Ca) a 2 ln/xr, (17) 

z i% = ^C F a s In fir, (18) 

where again all the derivatives are assumed to act on the right. As discussed in Sec. 3 the 
potential (in particular the matching coefficients given in Eqs. (7)-(14)) is ambiguous due to 
the freedom to perform time-independent unitary field redefinitions. This ambiguity affects 
also the normalization factor, Z X J 2 , so that terms in the potential can be traded for terms 
in Z l J 2 in a correlated way. Moreover, we note that, once the potential has been fixed, the 
normalization factor maintains at least a residual ambiguity, which amounts to iO, O being 
a Hermitian operator that commutes with the singlet Hamiltonian, p 2 jm + V} ' + . . .. 

Finally, let us note that the order l/m° and 1/m terms seem to be protected against hard 
contributions O(lnm). This has been explicitly proved for the l/m° term (the static poten- 
tial) at any finite order ina s [19]. 



5 Heavy quarkonium spectrum 



In this section we compute the heavy quarkonium spectrum up to, and including, the leading 

logs at 0{ma>l). Let us put ourselves in the situation Aq CD < mv 2 . The size of each term in 

the pNRQCD Lagrangian (2) is well-defined and can be evaluated as follows. The relative 

momentum p, and the inverse relative coordinate 1/r have a size of 0(ma s ). The time 

derivative has a size of 0(ma 2 ). US gluon fields, derivatives acting on it and the centre of 

mass momentum (in the rest frame, when entering in recoil corrections, due to the virtual 

emission of US gluons) have a size of the order of the next relevant scale (ma 2 or Aqcd)- 

Therefore, in order to obtain the leading logs at 0(ma s 5 ) in the spectrum, VJ ^ has to be 

computed at 0(a* In), U s (1) at O(afln), U s (2) at 0(a s 2 ln) and V s {3) at 0(a s ln). From the 

previous matching, we have obtained VJ°\ V} 1 ^ and VJ 2 ^ with the desired accuracy. Since 

terms of the type 0(a s In) can not exist at tree level, VJ 3 ^ does not need to be computed. The 

only term to be considered at 0(1 /m 3 ) needs to be a s -independent and it is the correction 

A 2 
to the kinetic energy, — - — -. 

Am 6 

After these considerations we can obtain the 0(ma^) leading-log correction to the heavy 
quarkonium spectrum. From the potential-like terms (4)- (6) we obtain the following cor- 
rection to the energy: 
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The lno; s appearing in Eq. (19) come from logs of the type In . Therefore, they can be 



mr 



deduced once the dependence on lnm is known. We have checked that our dependence on 
lnm coincides with the one obtained in Ref. [15]. The \x dependence of Eq. (19) cancels 
against contributions from US energies. It agrees with the \x dependence found in [14] for 
/ = 0. At the next-to-leading order in the multipole expansion, the contribution from these 



scales reads 



S Vb E n ^) 



3N, 



T F /rft(n,/|re^"- H °)r|n,/)(E a (t)0(t,O)^ j E 6 (O))(/i), (20) 



where H = \- V} ' and a is the UV cut-off of pNRQCD. Then, the total correction to 

2m 
the energy is 



SE n<u = 5 pot E n ^(fi) + 5 vs E ntl (fi) . 



(21) 



Different possibilities appear depending on the relative size of Aqcd with respect to the US 



scale ma 2 . If we consider that A QCD 



ma: 



the gluonic correlator in Eq. (20) cannot be 



computed using perturbation theory. Therefore, in a model independent approach, one can 
leave it as a free parameter and fix it with an experiment at some scale fi (since the running of 
Eq. (20) with /i is known, one can then obtain its value at another scale). Another possibility 
is to try to obtain it from lattice simulations (see for instance [20]) or by using some models 

[21]. 

If we consider that ma 2 ~^> Aq CD , Eq. (20) can be computed perturbatively. Since ma 2 is 
the next relevant scale, the effective role of Eq. (20) will be to replace jj, by ma 2 (up to finite 
pieces that we are systematically neglecting) in Eq. (19). Then Eq. (21) simplifies to 
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(22) 



Since in this situation one is assuming that 



Aqcd 
ma 2 



<^i 1, one can expand on this parameter. 



Therefore, non-perturbative corrections can be parameterized by local condensates. The lead- 
ing and next-to-leading non-perturbative corrections have been computed in the literature 

[22,23]. 



When working in the pole-mass scheme, it is expected that large corrections will appear from 
the finite pieces at NNNLO. This problem seems to be solved by using renormalon-based 
mass definitions or the like [24]. In this case we could expect our calculation to provide an 
estimate of the magnitude of the NNNLO corrections. 



6 Conclusions 



We have computed the matching between NRQCD and pNRQCD and the heavy quarkonium 
spectrum at the leading-log NNNLO. Within the pNRQCD effective field theory framework, 
our results almost trivially follow from existing calculations. The results we achieved are 
important at least for t-t production and T physics. In the first case the results we present 
are a step towards the goal of reaching a 100 MeV sensitivity on the top quark mass from 
the t-t cross-section near threshold to be measured at a future Next Linear Collider [25]. In 
the second case it will improve our knowledge on the b mass. 

Our calculation paves the way for a complete NNNLO analysis of heavy quarkonium. There- 
fore, let us briefly comment here on what extra calculations are still required in order to 
obtain the pNRQCD Lagrangian that is relevant at NNNLO and the spectrum at 0(ma^). 
The contribution from the US scales have been computed in the MS scheme at one loop 
[14]; also the matching between QCD and NRQCD is known at the one-loop level in the 
MS scheme [2,12]. Therefore, the missing calculations concern only the matching between 
NRQCD and pNRQCD. The static potential needs to be known at three loops, the 1/m 
corrections at two loops, the 1/m 2 at one loop and the 1/m 3 (if any) at tree level, all of them 
in the MS scheme so as to profit from the already known results. 
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